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We study the S = 1 square lattice Heisenberg antiferromagnet with spatially anisotropic nearest 
neighbor couplings Ji x , Ji y frustrated by a next-nearest neighbor coupling J2 numerically using the 
density-matrix renormalization group (DMRG) method and analytically employing the Schwinger- 
Boson mean-field theory (SBMFT). Up to relatively strong values of the anisotropy, within both 
methods we find quantum fluctuations to stabilize the Neel ordered state above the classically 
stable region. Whereas SBMFT suggests a fluctuation-induced first order transition between the 
Neel state and a stripe antiferromagnet for 1/3 < Ji x /Ji y < 1 and an intermediate paramagnetic 
region opening only for very strong anisotropy, the DMRG results clearly demonstrate that the 
two magnetically ordered phases are separated by a quantum disordered region for all values of the 
anisotropy with the remarkable implication that the quantum paramagnetic phase of the spatially 
isotropic J1-J2 model is continuously connected to the limit of decoupled Haldane spin chains. Our 
findings indicate that for S = 1 quantum fluctuations in strongly frustrated antiferromagnets are 
crucial and not correctly treated on the semiclassical level. 

PACS numbers: 75.50.Ee, 75.10.Jm, 75.30.Kz 



I. INTRODUCTION 

A striking commonality between the recently discov- 
ered iron pnictides superconductors'^^ and the high- 
T c cuprates is that in both cases superconductivity 
emerges on doping antiferromagnetic parent compounds. 
In trying to unravel the mechanisms of superconductiv- 
ity frustrated quantum antiferromagnets on the square 
lattice have been subject to intense research over the last 
decades with particular interest in the extreme quantum 
limit S = i relevant for the cuprates, whereas the inter- 
est in higher spin values increased tremendously with the 
discovery of the pnictides. 

Thereby, the Heisenberg model with antiferromagnetic 
exchange couplings Jj and J2 between nearest neighbor 
(NN) and next-nearest neighbor (NNN) has served as a 
prototype model for studying magnetic frustration. On 
a classical level, one finds Neel order to be stable for 
J2/J1 < 1/2 whereas for larger ratios the classical ground 
state is give by a stripe-antiferromagnet with ordering 
wave- vector (tt, 0). Not surprisingly, the J1-J2 model has 
been used to rationalize the (tt, 0) magnetism 6 - of the 
iron pnictide superconductors and to subsequently cal- 
culate the magnetic excitation spectra 7 -* 8 - where the in- 
corporation of a strong anisotropy between the NN cou- 
plings turned out to be necessary to reproduce the low 
energy spin- wave excitations. Recently, it has been sug- 
gested that the strong anisotropies in the magnetism 8 , 
but also in electronic properties 9 - of the pnictides origi- 
nate in the coupling to orbital degrees of freedom arising 
from an orbital degeneracy of an intermediate 5=1 
spin stated. Even if the charge degrees of freedom are 



not completely localized, the J\-Ji model can be used 
as the starting point for a symmetry-based analysis of 
magnetism and superconductivity in the iron-based su- 
perconductor o 11 ' 12 . 

The 1/S expansion serves as a natural starting point 
to investigate the stability of the classical orders against 
quantum fluctuations which are expected to induce a 
paramagnetic phase near J2/J1 = 1/2 where both or- 
ders compete. Indeed, on the level of lowest order lin- 
ear spin-wave theory one finds an intermediate param- 
agnetic phase for all spin values 1 ^. However, the in- 
corporation of spin-wave interactions within a modified 
spin-wave theory (MSWT)! 4 -^ or the SBMFTi 6 - drasti- 
cally changes this picture. Within both approaches one 
finds a dramatic stabilization of the classical orders by 
quantum fluctuation for the Neel order even up to val- 
ues considerably larger than the classical threshold value 
J2/J1 = 1/2. This has been interpreted as an order-out- 
of-disorder phenomenon giving rise to a fluctuation in- 
duced first order transition between the two orders. This 
picture has been confirmed recently by a functional one- 
loop renormalization group analysis 1 - 7 - where it was shown 
that the Neel phase becomes unstable towards a fluctua- 
tion induced first order transition for S > 0.68. 

Surely, it is questionable if the semiclassical treatment 
can correctly account for the quantum fluctuations in the 
regime of strong frustration and small spins. In fact, in 
the case S = \ various numerical studies including exact 
diagonalizatio n 18 i 19 ' 20 i 21 , variational Monte Carlo 22 i 23 se- 
ries expansio n 24 ' 25 ! 26 ' 27 ! 28 as well as the coupled cluster 
approach 29 give the consistent picture that in the regime 
0.4 < J2/J1 < 0.6 no magnetic order is present clearly 
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indicating that the aforementioned semiclassical treat- 
ments overestimate the stability of the ordered states. 
Another key observation of the series expansion and in 
particular of the unbiased exact diagonalization studies 
is that in the paramagnetic phase the lattice symmetry is 
spontaneously broken due to the formation of columnar 
valence bond solid order. Such states have been pre- 
dicted befor o 30 ' 31 on the basis of a large N treatment for 
non frustrated systems in the absence of long-range an- 
tiferromagnetic order. These depend however crucially 
on the value of the spin, revealing the subtle workings 
of the spin Berry phases (see e.g. [HJ and references 
therein). For S — \ these act to give a special stability 
to 'valence bond' pair singlets involving nearest-neighbor 
spins, and the quantum disordered phases are actually 
valence bond crystals breaking the translational symme- 
try of the square lattice further to fourfold degenerate 
spin-Peierls states. But for S = 1 the Berry phases act 
in favor of the formation of 'chain singlets ' 30 ' 31 that can 
stack either along the x or y direction of the square lattice 
yielding a twofold degenerate ground state. 




FIG. 1: The model studied in this paper consists of a S = 1 
Heisenberg system with spatially anisotropic nearest-neighbor 
couplings Ji x , Ji y and isotropic next-nearest-neighbor cou- 
plings J2. Pending the balance of these couplings one finds 
either a simple staggered- (a) or 'stripe' (b) antiferromag- 
netic order. In the limit of isolated chains a stack of isolated 
Haldane spin chains is formed (c) and based on our DMRG 
calculations we could conclude that these survive all the way 
to the isotropic limit in the vicinity of the point of maximal 
frustration (J2/J1 = 0.5). 

In this paper we investigate the S = 1 version of the 
spatially anisotropic J\-J% model, given by the Hamilto- 
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where the first two sums run over NN spins with exchange 
couplings Ji x < Ji y in x and y directions, respectively, 
and the third sum runs over all NNN pairs with exchange 
couplings J2 as illustrated in Fig [T] We note that both 
the spatially isotropic J\-Ji model as well as the decou- 



pled Haldane spin-chain limit appear as special cases of 
Hamiltonian |T]). 

Unlike for S — | this model has been hardly explored 
for S — 1 despite the potential relevance for the iron 
pnictides, and to best of our knowledge exact diagonal- 
ization studies are not available. In a very recent cou- 
pled cluster treatment^ 3 - it has been found that at the 
isotropic point J\ x = J\ y = J\ the Neel and stripe or- 
dered phases are separated by a first order transition at 
J2/J1 « 0.55 slightly smaller than the MSWTi 4 .^ and 
SBMFT— results and that the transition remains first 
order up to an anisotropy J\ x j J\ v ps 0.66. For stronger 
anisotropies the authors find continuous transitions very 
close to the classical transition line J2 = Ji x /2, although 
they were not able to resolve an intermediate paramag- 
netic region within the numerical resolution. On con- 
trary, a two step DMRG study 3 ^ at the particular point 
Jix/ Jiy — 0.2 indicated a much wider non-magnetic re- 
gion with a spin-gap reaching its maximum A pa 0.39 J± y 
close to the maximally frustrated point J 2 / ' Ji y = 0.1, 
which is only slightly below the gap Ah ~ 0.41J [35j of 
an isolated Haldane chai n 36 ' 37 . 

Starting from the Haldane chain limit ( J\ x = J2 = 0) it 
has been predicted by Monte Carlo simulations^ as well 
as analytically 3 - that an infinitesimal coupling J\ x leads 
to the destruction of the topological string order— of the 
Haldane chain. However, due to the protection by the 
finite energy gap of the spin-1 Haldane chain all other 
ground-state properties as well as the thermodynamics 
are only minimally affected by a small interchain coupling 
J\ x and a finite, albeit smal l 41 ' 42 ' 43 ' 44 coupling is neces- 
sary to establish magnetic long range order. Whereas a 
similar reasoning should hold for a vanishing NN cou- 
pling perpendicular to the chains (J\ x — 0) and small 
diagonal coupling J2 the Haldane chain phase seems to 
be considerably more stable against the simultaneous in- 
crease of the two mutually frustrating couplings J\ x and 
J2 as indicated by two step DMRG results showing an 
almost negligible reduction of the spin gap for moder- 
ate interchain couplings 34 . This suggests that the para- 
magnetic phase of the frustrated two dimensional spin 
model can be viewed as a continuation of the Haldane 
spin chain. Recently, based on a theoretical method 
with a continuous deformation of the J\-J% model, the 
ground state at the special isotropic case J2 = Ji/2 has 
been conjectured 4 ^ to be a two-fold degenerate valence 
bond solid state along either the horizontal or vertical 
direction of the square lattice. The main result of this 
paper is that by employing DMRG and by studying sys- 
tematically how matters evolve as function of increas- 
ing anisotropy we arrive at solid evidence for the smooth 
continuation between the Haldane chain phase and an 
isotropic disordered phase near J 2 = Ji/2. 

Since for S — 1 exact diagonalization is restricted to 
very small system sizes and since the quantum Monte 
Carlo method suffers from the infamous sign problem 
in the presence of frustration we employ the DMRG 
method^ 6 - to map out the phase diagram of the spatially 
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anisotropic J\ x —J\y—Ji model (QJ over the whole param- 
eter range including the decoupled Haldane spin chains 
and the isotropic J1-J2 model as limiting cases. The 
DMRG has the merits of not being biased as analyti- 
cal treatments based on 1/S or 1/N expansions, being 
capable of dealing with frustrated systems, and allow- 
ing to investigate much bigger system than possible with 
exact diagonalization. Moreover, it is capable of repro- 
ducing the spin gap of the decoupled Haldane chain limit 
with high accuracy^. The resulting phase diagram is 
contrasted by analytical results we obtain within a gen- 
eralization of the SBMFT to the anisotropic system. 

The remainder of the paper is organized as follows. In 
Sec.|ll]wc outline the SBMFT and calculate the resulting 
phase diagram which is used for comparison with the nu- 
merical results. In Sec. lIIII we use the DMRG to construct 
the phase diagram by a careful finite-size analysis of the 
magnetic structure factor at the ordering wave vectors 
of the two magnetically ordered phases as well as of the 
spin-gap in the non magnetic region as a function of the 
NNN coupling J2 for various ratios of the NN couplings 
J\ x and Ji y . Finally, in Sec. [IV] the results obtained 
within the two methods are compared and the shortcom- 
ings of the semiclassical approach are highlighted. Fur- 
ther, the potential relevance of our findings for the iron 
pnictides is discussed. 



where the sum runs over all bonds and we have in- 



troduced the bond operators F i 



b\ _ v . A mean-field decoupling is then performed with 

respect to the order parameters — (FjA/2 and gij = 

{G[-)/2. For the Neel ordered phase we have to introduce 
fields g x , g y for the non- frustrated nearest-neigbor bonds 
and fd for frustrated NNN bonds whereas the (0, ir)- 
phase is characterized by order-parameter fields g y and 
gd for the non-frustrated bonds along the y and diag- 
onal directions and f x for the frustrated x-bonds. The 
local constraints are replaced by a global one and treated 
with a Lagrange multiplier A. The resulting mean-field 
Hamiltonian, which in momentum space is given by 



6 i,A" and G ij 



t _ 



{K + A)(^ T & q1 + bl ql b- ql ) 



(4) 



is easily diagonalized by a Bogolioubov transformation 
yielding the dispersion 



u q = [{h q + A) 2 A 2 // 2 



(5) 



with 



II. SBMFT 

In this section we generalize the SBMF calculation for 
the isotropic J\- J2-modeli£ to the case of anisotropic 
nearest-neigbor exchange couplings. For abbreviation we 
define J x :— J\ x , J y := J\ y , and Jd '■= Ji- Assuming 
Jx < Jy, the classical ground states are given by a Neel 
ordered state with an ordering wavevector Q — (tt, tt) for 
Jd/Jx < 1/2 and by a columnar antiferromagnetic state 
with Q = (0, tt) for Jd/ Jx > 1/2. Following the previ- 
ous calculations, we perform a spin rotation Sf = <TiSf , 
Sf = aiS?, Sf = Sf, where a, = exp(iQn) = ±1 
with Q the ordering wavevectors of the two classical or- 
ders, and represent the rotated spin operators in terms 
of Schwinger bosons bj j, bi_i as 



(2) 



Here a = (a x , a v , a z ) with a a the standard Pauli ma- 
trices. The constraint b\ v bi^ = 2S ensures that {S 2 } = 
S(S + 1). The Hamiltonian ([1]) can then be rewritten in 
the compact form 



n = 



! ±±±(F} ]FlJ -2S 2 ) 



+ - 



-(GlG. 



ij l 3 



2S 2 ) 



(3) 



h q = &f d J d cos q x cos q y , (6a) 
A q = 2(g x J x cos q x + g y J y cos q y ) (6b) 

in the Neel phase and likewise in the (0, tt) phase 



h q = 2f x J x cosq x , (7a) 
Aq = 2g y J y cosq y + Ag d JdCOsq x cosq y . (7b) 

The Lagrange multiplier in the ordered phases is de- 
termined by the requirement that ojq = for the cor- 
responding ordering wavevectors yielding A = 2(g x J x + 
g y J y )~Af d Jd in the Neel and A = 2{g y J y - f x J x ) + Ag d Jd 
in the (0, 7r)-phase. The reduced moment S* in the or- 
dered phases is determined by 



s * = s _l( [h + x 

2 \Jq " g 



(8) 



Naturally, one might expect that quantum fluctuations 
tend to destabilize the classical orders introducing an in- 
termediate paramagnetic phase. In this case the second 
order transitions between the two different magnetic or- 
ders are determined by the lines where the corresponding 
magnetizations go to zero, S* — * 0. However, the anal- 
ysis of the isotropic modet^ (J x — J y ) shows that the 
quantum fluctuations can lead to a significant stabiliza- 
tion of the Neel order leading to a region where ^ > 
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and S( Q n j > 0. In this region which is expected to per- 
sist for not too strong anisotropy, a discontinuous first 
order transition between the two different magnetic or- 
ders is likely and can be estimated from a comparison of 
the ground-state energies which immediately follow from 
Eq. © as 



E, 



E, 



(tt.tt) = -Jx&gl - S 2 ) - J v (2g y - S) 

+2J d (2f 2 - S 2 ) (9a) 
'(0,-ir) = Jx(2f 2 - S 2 ) - J y (2g 2 - S 2 ) 

-2J d (2g 2 - S 2 ). (9b) 



The order parameter fields entering the mean-field 
Hamiltonian have to be determined self consistently. 
Using the above Bogoliubov transformation we obtain for 
the Neel ordered phase 



S 



A 



o — COSq x , 

q 2UJ Q 

S + / —2-cosq y , 

Jq ^q 



fd = s* 



(10a) 
(10b) 

j- cos q x cos q v , (10c) 

lLO„ 



h q + X 



whereas in the (0, w) phase the self-consistency equations 
for the order-parameter fields are given by 



stripe order 




Neel order 
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( J 1y" J 1x) /J 1y 

FIG. 2: Phase diagram for S — 1 as a function of and 
anisotropy a — (Ji y — Ji x )/ Ji y > between the nearest neig- 
bor exchange couplings and relative strength of the next near- 
est neighbor coupling J2 / Ji y obtained within SBMFT. Up to 
an anisotropy a « 0.66 the (n, n) and (0, tt) antiferromagnetic 
orders are separated by a first order transition (solid line) . At 
the tricritical point the first-order line splits into two second- 
order lines (dashed and dotted) separating the two magnetic 
phases from a gapped non-magnetic phase. 



fx = S*+ [ ^L±±cosq x , (11a) 

Jq Zu > q 

g v = S * + J ^j~ cosq y> ( llb ) 

g d = S* + [ ^-coaq x cosq y . (11c) 

Jq l ^q 

The resulting phase diagram for J\ x < J\ v is shown in 
Fig. [2] as a function of the anisotropy a = {J\ v — Jix)/ Jiy 
between the NN exchange couplings and of the relative 
strength (3 = J2/ 1 J\ y of the NNN coupling. In agreement 
with earlier studies of the isotropic limit J\ x = J\ y = J\ 
within SBMFT and MSWT we find a dramatic stabi- 
lization of the Neel order above the classical value and 
a large region 0.51 < j3 < 0.68 where the two compet- 
ing orders are potentially stable, indicated by SL ^ > 
and S* w <. > 0. The crossing of the self-consistently de- 
termined energies of the two states suggests a first or- 
der transition at P ss 0.64 considerably larger than the 
classical value 1/2 and about 10 percent bigger than the 
coupled cluster result 33 . Although the region of coex- 
istence is considerably narrowed by a small anisotropy 
we find it to persist up to a 0.66 where the first or- 
der line terminates. The existence of such a tricritical 
point was also suggested by the coupled cluster analysis 
although located at a much smaller anisotropy a « 0.34. 



For anisotropics a > 0.66 we find an intermediate non- 
magnetic region (S 1 ^ ^ = S* Q ^ = 0) separating the two 
ordered phases. This region is found to be very narrow 
and close to the classical phase boundary (3 — (1 — a)/2. 

Although the existence of a tricritical point and of a 
very narrow paramagnetic strip terminating at the Hal- 
dane chain limit J\ x = J2 — is in qualitative agreement 
with the coupled cluster results^ 3 , the small width of the 
paramagnetic region is in disagreement with other avail- 
able numerical results. For J\ x j J\ y = 0.2 (a — 0.8) a 
two step DMRG calculation^ shows a much wider non 
magnetic region centered around the classical transition 
point J-z/Jiy = 0.1. Interestingly, the spin gap at this 
maximally frustrated point is almost identical to that of 
an isolated Haldane chain suggesting that even for rela- 
tively strong interchain couplings one dimensional Hal- 
dane chain physics is still important. This is certainly 
missed by the SBMF calculation which treats the spin as 
a continuous variable and does not distinguish between 
integer and half-integer spins crucial for the existence 
of the Haldane spin gap. Moreover, it has been estab- 
lished within quantum Monte Carlo calculation a 41 ' 42 i 43 
that for J 2 = the transition between the Neel ordered 
phase and the gapped non-magnetic phase is located at 
Jix/ Jiy = 0.044, again indicating that the paramagnetic 
region close to the Haldane chain limit is considerably 
wider than suggested by both the SBMFT and the cou- 
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pled cluster results^ 3 .. Since for J2 — the system is not 
frustrated quantum Monte Carlo can be considered exact 
in this regime. 

We have also calculated the transition out of the (0, if) 
state in linear spin-wave theory, for comparison to the 
SBMFT calculation and to preliminary neutron scatter- 
ing results on the pnictides that found (0, tt) order but 
with a small moment. This method computes the re- 
duction of the classical antiferromagnetic moment due to 
zero-point excitations of spin waves (which captures the 
1/5 correction to the classical moment in a large- S ex- 
pansion) . The transition line out of the ordered phase is 
estimated as the point where the correction is as large as 
the original moment. 

We assume three antiferromagnetic couplings: J\ x , 
Ji y , and J2 with J\ y > J\ x and that we are in a (0, it) 
ordered phase, which requires (as is evident from the for- 
mula below) that J\ x < U^. The dispersion relation 
for spin- wave excitations was previously obtained for this 
anisotropic J\-J% model in Ref. [l^- The integral for the 
correction to the classical moment around the (0, tt) case 
is, in units of the Bohr magneton and with lattice spacing 
a = 1, 



Am 



wit : ] 



d 2 k 



[0,2^ (2tt) 2 cos 2 (k y ) I f(k x y 



- 1 



2J 2 + Jly - Ji x [l - cos(fc x )] 
2J 2 cos(k x ) + Jiy 



(12) 



(13) 



The transition is found numerically to lie very close 
to the classical transition line J\ x — 2J 2 : the normal- 
ized difference (Jf — J\ x /2)/ Ji y , where J$ is the criti- 
cal coupling where the correction is equal to the original 
moment, is always less than 2% for < J\ x < 0.99Ji a . 
Significant reduction of the ordered moment also occurs 
only near the classical transition line. The primary dif- 
ference from the SBMFT calculation is that the spin- 
wave calculation always predicts a paramagnetic phase 
between stripe and Neel order. The width of this param- 
agnetic phase is much smaller in cither approach than in 
the DMRG calculation of the following section, because 
these analytical approaches do not capture the strong 
quantum fluctuations that favor the Haldane phase. 



indicates that the SBMFT tends to overestimate the sta- 
bility of the magnetically ordered phases, surely close to 
the Haldane chain limit but presumably also for larger 
values of J\ x and J 2. In the following we shall refine the 
boundary region by using DMRG method. 

In the following DMRG calculation, we will set J% y = 
Ji as energy unit, and a periodic boundary condition 
(PBC) is used and in each DMRG block up to m = 3200 
states are to be kept with the truncation error in the 
order of or less than 10 -5 . 



A. Isotropic case with J\ x — J\ y 



Ji 



Let us first consider the isotropic case J\ x = J\ v — 
Ji where the SBMFT suggests the strongest tendency 
towards a first order transition although the transition 
point J2/J1 ~ 0.64 obtained from a comparison of the 
energy minima seems to be suspiciously high compared to 
the classical transition point J2/J1 = 0.5. Fig. [3] shows 
the ground state energy per site calculated by DMRG 
with the sample size varying from N = 4 x 4 (16 sites) 
up to N — 8 x 8 (64 sites). The ground state energy 
reaches the maximum with J 2 between 0.54Ji and 0.58 Ji, 
which becomes sharper with the increase of sample size, 
indicating a region with possible phase transitions below 
the first order transition point obtained in SBMFT but 
still notable above the classical transition. 




0.4 0.6 
J2/J1 



FIG. 3: (color online) The ground-state energy per site at 
different system size A = 4 x 4 (blue triangle), 6x6 (red 
circle) and 8x8 (black square), for the isotropic case with 

Jim — Jly — ^1 • 



III. DMRG 



We examine such a region by calculating the magnetic 
structure factor S z (q), defined by 



The above SBMFT has clearly shown an interesting 
narrow boundary region between the Neel and stripe or- 
dered phases in the phase diagram of Fig. [21 where the 
quantum fluctuations are expected to become very im- 
portant. Most interestingly, it suggests an increasing ten- 
dency towards a fluctuation induced first order transition 
on approaching the isotropic point J\ x = J\ y . However, 
the comparison with previous numerical result o 34 i 41 ' 42 i 43 



-iq(rj-rj) 



(SfSj) 



As expected, we find that S z (q) shows a dramatic change 
from peaking at Q = (tt, tt) (the Neel order) to Q — (0, ir) 
(the stripe order) with increasing J2/J1. Fig. |U (a) and 
(b) illustrate S z (q)/N vs. J2/J1 for system sizes N — 
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FIG. 4: (color online) The evolution of the peak values of the 
structure factors, S z (tv,tt)/N (a), S z (0,n)/N (b), and of the 
spin-1 gap (c), at three different size JV = 4x4, 6x6 and 
8 x 8, as well as their thermodynamic limit extrapolations in 
the isotropic case Ji x = J\ v = Ji = 1. The finite-size scaling 
for the spin-1 gap is shown in (d). 

4 x 4, 6 x 6 and 8 x 8, as well as the thermodynamic limit 
values obtained by a finite size scaling. Here we have 
used the quadratic function f(x) = A + Bx + Cx 2 with 
x = j7 to perform the finite-size scaling. 

The Neel order transition point is found at J2/J1 ~ 
0.525 in Fig. 0Ja), which is clearly distinct from the 
stripe order tranistion point at J2/J1 ~ 0.55 in Fig. BJb) 
indicating that the two magnetically ordered phases are 
separated by an intermediate non-magnetic region. 

To independently verify the above results, we also cal- 
culate the spin-1 gap AE(S = 1) = Ei(S = 1) - E 
presented in Fig. @fc) and (d). Fig. |3Jd) illustrates 
the finite-size scaling for the spin gap AE(S = 1) at 
different values of J2 using a scaling function f(x). In 
Fig. BJc), the evolution of the spin-1 gap as a function of 
Ji is given. In the thermodynamic limit, the transition 
points determined by the spin-1 gap are at J2 ~ 0.525 
and J2 ~ 0.555, respectively, which are very close to 
the previous results determined by the structure factor. 
Therefore, for the present 5 = 1 J1-J2 model in the 
isotropic limit, our numerical approach has established 
an intermediate spin disordered region with a finite spin 
gap which separates the two ordered magnetic phases. 



B. Anisotropic case with J\ x < J% y = Ji 

Now we consider how the spin disordered phase evolves 
with the increase of anisotropy at J\ x < J\ y = J\. First 
we consider the case at J\ x = 0.5Ji, and the results are 
presented in Fig. [5J By using the same finite size scal- 



FIG. 5: (color online) The peak values, (a) Sz(tt,tv)/N, 
(b) S z (0, n)/N, and (c) the spin-1 gap, versus J%/Ji v at 
Jix/Jiy = 0.5 for different sizes. The finite-size scaling for 
the spin-1 gap is shown in (d). 

ing procedure, we find that the spin disordered phase is 
bound by a lower transition point J2/J1 ~ 0.24 and an 
upper transition point J2/J1 ~ 0.28 based on the struc- 
ture factor calculation. Again the spin gap calculation 
gives rise to a consistent spin disordered regime between 
J2/J1 « 0.23 and J 2 /J x » 0.285 as shown in Fig. [51(c). 

In the same regime, with a fixed J2/J1 — 0.25, we 
have further studied the phase boundaries by varying 
Jix/Ji- As shown in Fig. [6l the lower transition point 
is obtained at J\ x /J\ — 0.45 and the upper transition 
point at Jix/Ji ~ 0.52, while the finite size scaling 
for the spin-1 gap results in a similar region between 
J lx /J x = 0.44 - 0.54. 

In the extreme case at J\ x = and J2 = with 
Jiy = J\ 1 the system simply reduces to an array of de- 
coupled S = 1 spin chains with a finite Haldane gap. 
Fig. [7] shows how the ground state continuously evolves 
from that of the well-known decoupled spin chains to the 
anisotropic 2D case by turning on J\ x / J\, which remains 
disordered until the Neel order sets in at J\ x j J\ ~ 0.05 
with S z (tt,tt)/N becoming finite in the thermodynamic 
limit. Indeed the corresponding spin-1 gap continuously 
decreases with the turning on of a finite J\ x , but only 
vanishes around J\ x j J\ ~ 0.06 as shown in Fig. [7] (c) 
and (d). 

Now we turn on J2- At J\ x = 0, we find that while the 
calculated structure factor at (n, n) continuously reduces 
as the sample size increases from TV = 4 x 4, 6x6, to 8x8, 
and is extrapolated to zero by finite size scaling, a finite 
stripe order S z (0,n)/N will emerge at J2I J\y = 0.025 
in the thermodynamic limit, which is further supported 
by vanishing spin-1 excitation gap at the same point as 
shown in Fig. [5) It is noted that the best finite scaling 
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FIG. 6: (color online) (a) &(tt, 7t)/jV ; (b) S z (0,ir)/N, and (c) 
the spin-1 gap, versus Ji x /Ji y at J2/J\ y = 0.25 for different 
sample sizes. In (d), the finite-size scaling for the spin-1 gap 
is given. 



FIG. 8: (color online) At J\ x = 0, the structure factors 
Sz(tt,tv)/N (a) and S z (0,tt)/N (b) as well as the spin-1 gap 
(c) versus J2 / Ji y are shown at different sizes, including their 
thermodynamic extrapolations. 
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FIG. 7: (color online) (a) S z (tt,it)/N, (b) S z (0, n)/N, and (c) 
the spin-1 gap, versus Ji x /Ji y at J2 = at different sizes. In 
(d), the finite-size scaling for the spin-1 gap is given. 



for the structure factor here is obtained with using a f(x) 
with x = 1/VN = 1/Ny, for a square lattice N x = N y , 
instead of x = 1/N used previously. The justification 
of such a finite-size scaling for the spin structure factors 
at Ji x = and small J2I J\y is given in Fig. [5] and its 
caption. 
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FIG. 9: (color online) (a) The size dependence of the struc- 
ture factor [for a given N y — 6 (6-leg systems)] at J\ x = and 
small J2/ 'Ji y which quickly saturates at N x > N y , indicating 
that x = 1/N y in the scaling function f(x) for the structure 
factor is more appropriate in this extreme limit. The corre- 
sponding finite-size scaling of the structure factors for square 
samples are illustrated in (b) and (c), respectively. For com- 
parsion, the spin-1 gap at different J 2 is still well scaled by 
x — 1/N in the scaing function as shown in (d). 
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C. Phase diagram 

The resulting phase diagram obtained by the DMRG 
calculations with careful finite size scalings of the mag- 
netic structure factor as well as the spin gap is shown 
in Fig. |TD] as a function of the anisotropy a = {J\ y — 
Jix)/J\y between the NN couplings and the relative 
strength (3 — J-il J\y of the NNN superexchange. 

In contrast to the SBMFT which in the regime of small 
anisotropy a suggests a first order transition between the 
Neel and stripe ordered phases (see Fig. [2]), the DMRG 
results clearly indicate a paramagnetic strip separating 
the two magnetically ordered phases for all values of the 
anisotropy a including the isotropic point J\ x — J\ y (a — 
0). With increasing a the width of this region is found 
to slightly increase. Close to the Haldane chain limit the 
SBMFT predicts an intermediate non magnetic phase, 
though the width of this region is found to be much larger 
in the DMRG calculation in quantitative agreement with 
a previous two step DMRG analysis at fixed a = 0.8 
(Ref.IH). 

For J2 = we can compare our results to recent QMC 
simulation a 41 ! 42 ' 43 , which in this regime can be consid- 
ered exact since the system is not frustrated and QMC 
is therefore free of any sign problem. Within DMRG we 
find the transition between the gapped Haldane phase 
and the Neel ordered phase to occur at J\ x j J\ y ~ 0.05 
(see Fig. ITUl) in good agreement with QMC transition 
point Jix/Jiy — 0.044. Moreover, in the Haldane chain 
limit J\ x — J2 = we obtain a spin gap A w 0.4 (see 
Fig. Efd)) very close to the exact value A# = 0.41 for 
the Haldane spin chain™^, again demonstrating that the 
finite size scaling is well converged. 

Very interestingly, starting from the Haldane chain 
limit we find the maximum spin gap in the paramagnetic 
phase to decrease only very slowly with increasing cou- 
plings Ji x and J2. Up to J\ x j J\ y = 0.5, the gap is almost 
identical to the Haldane spin gap (see Fig. EJd)) indicat- 
ing that the Haldane spin chain is very robust against the 
simultaneous increase of the mutually frustrating cou- 
plings J\ x and J2. Remarkably, even the paramagnetic 
phase of the isotropic J1-J2 model with a maximum spin 
gap still being about half of the Haldane gap is continu- 
ously connected to the Haldane spin-chain limit. 

At the isotropic point J\ x = J\ y (a = 0) we 
find an intermediate paramagnetic phase for 0.525 < 
J2/J1 i$ 0.555. This region is considerably smaller 
than in the S = h case, where a paramagnetic regime 
0.4 < J2/J1 < 0.6 presumably with columnar valence 
bond order has been established by various numeri- 
cal m aho&&&&2L2L®&&&£L22.. Interestingly, the 
stabilization of the Neel phase above the classical transi- 
tion point is in agreement with the SBMFT. 

Despite the apparent failure of the semiclassical 
SBMFT in dealing with the strong quantum fluctuations 
in the boundary region between the two magnetically or- 
dered phases it is interesting to note that the param- 
agnetic region covers the phase boundaries obtained by 
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FIG. 10: (color online) Phase diagram for the anisotropic J%- 
J2 model determined by DMRG. The Neel order and stripe 
order phases are separated by a paramagnetic regime with 
the boundaries denoted by the red line with full circles for 
the upper and the blank line with full squares for the lower 
transition points, respectively. In the middle of the param- 
agnetic region lies in a dotted line with crosses which marks 
the maximal spin-1 gap AiJ max . Note that the symbols of cir- 
cle, square, and triangular denote the phase transition points 
determined by DMRG in the previous figures. 



SBMFT for the most of J\ x < J\ y = J\ region except for 
the part close to the isotropic limit and that the points 
where the spin gap is maximum (indicated by crosses and 
dotted in Fig. [IT))) are almost on top of both the first or- 
der transition line as well as of the two narrow second 
order lines obtained within SBMFT. 



IV. DISCUSSION 

In summary, we have studied the frustrated spin-1 
Heisenberg J\ x -J\ y -J2 model on a square lattice numeri- 
cally using the DMRG method and analytically employ- 
ing the Schwinger-Boson mean-field theory (SBMFT). 
Interestingly, this model contains both the isotropic J\-Ji 
model as well as decoupled Haldane spin chains as lim- 
iting cases. Moreover, it has attracted a lot of attention 
recently since it has been motivated as an effective model 
to describe the (0, 7r)-magnetism and the low-energy spin- 
wave excitations of the celebrated novel iron pnictide su- 
perconductors. Furthermore it has been suggested that 
the drastic reduction of magnetic moments to a value of 
0.4 hb is caused by strong quantum fluctuations in the 
vicinity of a continuous phase transition. However, to the 
present day the phase diagram has not been determined 
yet. 

Within both the SBMFT and the DMRG we find that 
the Neel phase is stabilized considerably by quantum fluc- 
tuations above the classically stable region up to a rela- 
tively strong anisotropy between the NN couplings. How- 
ever, in all other regards the phase diagrams obtained by 
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the two methods clearly disagree, indicating the impor- 
tance of the strong quantum fluctuations. 

The SBMFT suggest a fluctuation induced first or- 
der transition between the Neel and the stripe antifer- 
romagnets terminating at a tricritical point and splitting 
into two second order lines separated by an intermedi- 
ate paramagnetic region only for large anisotropies. Al- 
though the existence of a tricritical point and a hardly 
seizable paramagnetic strip terminating at the decoupled 
Haldane spin chain point {J\ x — J2 — 0) are consistent 
with a recent coupled cluster calculation^ 3 ., the SBMFT 
definitely falls short on approaching the Haldane chain 
limit. This becomes clear by a comparison with quantum 
Monte Carlo simulation o 41 ' 42 ' 43 showing a much larger re- 
gion of stability of the Haldane chain phase against inter- 
chain couplings J\ x . Also the previous two-step DMRG 
calculation 34 at J\ x / J\ y — 0.2 indicates a paramagnetic 
region being an order of magnitude wider than found 
within the SBMFT. 

On the contrary, the phase diagram obtained within 
DMRG by a careful finite size scaling of the magnetic 
structure factor and the spin gap is consistent with both 
the QMC and the previous DMRG results. Furthermore 
the spin gap in the decoupled chain limit agrees well 
with the exact value for the isolated Haldane chain. The 
width of the paramagnetic region slightly decreases on 
approaching the isotropic point but remains finite over 
the entire parameter range. This has the remarkable im- 
plication that the paramagnetic phase of the frustrated 
two dimensional spin model, including the isotropic J\- 
J2 limit, is continuously connected to the Haldane spin- 
chain phase. In other words, the paramagnetic phase can 
be viewed as a continuation of the Haldane spin chain 
phase, with the caveat that the topological string order 
has to disappear at any finite interchain couplin g 38 i 39 . 

The reason for the failure of the semiclassical SBMFT 
in dealing correctly with the strong quantum fluctuations 
in the boundary region between the two classical orders is 
easy to understand. The SBMFT deals with the spin as 
continuous variable while it is blind for the Berry phase 
effects that distinguish between half-integer and integer 
spin values crucial for the existence of the Haldane gap as 
well the valence bond crystals. This physics is definitely 
missed by the semiclassical treatment. Our DMRG re- 
sults demonstrates that the Haldane chain phase is very 
robust against the simultaneous increase of the couplings 
J\ x and J2 along the strongly frustrated boundary re- 
gion, indicated by the minute reduction of the spin gap 
compared to the isolated Haldane chain limit. 

What does the study of this spin system teach us re- 
garding the superconductivity in the iron pnictides? One 
could speculate that the basic physics is similar as in the 
cuprates. Although the spins are larger in the pnictides, 
the geometrical frustration renders the spin system to 



be on the verge of undergoing a quantum phase transi- 
tion in a quantum disordered state. This quantum spin 
physics then sets the conditions for the emergence of su- 
perconductivity in the doped systems. But in this regard 
the size of the microscopic spin does matter more than 
one intuitively anticipates. For S = 1/2 the well estab- 
lished fact that the Berry phases conspire to turn the 
quantum disordered states of the spin-only systems into 
valence bond solids gives a rationale to take Anderson's 
resonating valence bond (RVB) idea for the origin of high 
Tc superconductivity quite seriously. The valence bonds 
are protected by the spin gap and the effect of doping 
could well be to just turn the valence bond solids into 
translational quantum liquids that transport two units 
of electrical charge. Focussing on the Berry phases in 
the pnictides one has only the options of an 'intermedi- 
ate' crystal field S — I stated or the high spin S = 2 
for the microscopic spin. In the former case the ground 
state has a twofold degeneracy and the 'building blocks' 
are no longer pair singlets but instead the 'chain sin- 
glets'. Thinking along the RVB lines, what to expect 
when such a system is doped? The dopants will increase 
the quantum fluctuations but chains are not pairs. One 
anticipates that doping might drive the system into the 
non-magnetic 'Haldane chain phase' breaking the two- 
fold rotational symmetry of the lattice. One notices that 
something of the kind is found in the phase diagram of 
the pnictides: the structural transition to the orthorom- 
bic phase persists to much higher dopings than the stripe 
antiferromagnetism 47 . At first view this seems rather de- 
tached from the RVB idea. But now one has to realize 
that the charge carriers are themselves spinfull, carrying 
by default a half-integer spin. Taking for instance the 
intermediate S = I background, the holes would carry 
S = 1/2, and t-J models corresponding with a mix of 
5=1 and S = 1/2 states have been studied in the 
past^. One anticipates that in an incompressible 'chain 
like' S = 1 background the S =1/2 carriers might again 
be 'glued by Berry forces' into valence bond pairs, re- 
establishing a connection with the RVB mechanism of 
the cuprates. A similar idea has been previously pro- 
posed from a different approach 4 ^. 

In conclusion, it is quite questionable that any of these 
considerations have a bearing on pnictide superconduc- 
tivity but they do make the case that there is still much 
interesting physics to be explored dealing with systems 
characterized by a larger microscopic spin. 
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